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Introduction The upper atmospheres of short periodic ex- 
oplanets are affected strongly by the stellar X-ray an d EUV 
(XUV ) radiation and plasma environment (Schnei der et al. I 
Il998l). Therefore, these b o dies should experience high ther- 



mal ( Lammer et al . ' '2003t_ Lecavelie r des Etangs et al.1 2004 



Yelle 



2004t 



Vid al Madiar et al. I 



2003; Baraffeetal. 



2004 



Tian et a l. 2005; Erkaev et al. 2006; Penz et al. 2006) and non- 
thermal (Erkaev et al. 2005; Khodachenko et al. 2006) mass 
loss. Khodachenko et al (2006) studied the minimum and 
maximum possible atmospheric erosion of the "Hot Jupiter" 
HD209458b due to Coronal Mass Ejections (CMEs) and found 
that this exoplanet, which orbits a solar-like star at 0.046 AU 
could have been eroded to its core-mass if its atmosphere were 
not protected by a strong magnetic field. Because the mass loss 
depends on the strength of planetary magnetic fields, which are 
estimated for slow rotating tid ally locked "Hot Jupite rs" to be 
in a range of 0.005 - 0. 1 Mjup dGriefimeier et al.ll2004l) . weakly 
magnetized short periodic gas giants can experience large non- 
thermal mass loss rates during their whole life time. The aim 
of this paper is to determine at which distances exoplanets can 
maintain their initial mass, and where close orbit gas giants can 
experience huge mass loss rates that may influence the short pe- 
riod CoRoT-planet population. 

1 . Exosphere-formation and the onset of planetary 
mass loss 

The overall evolution of a planet can be split into (1), the for- 
mation period during which the planet is formed and its mass 
grows, and (2), the mature period with the planet detached from 
gas reservoirs and exposed to loss processes driven by its host 
star. During the formation period the planet is in contact with a 
mass reservoir having certain pressure: the protoplanetary neb- 
ula. During an isolated period planets are surrounded by vac- 
uum with no continuum pressure and an exosphere as the inter- 
face. The overall mass history of a planet is thus determined by 
pressure gradients as long as there is a nebula pressure and by 
exosphere properties once there is no nebula pressure. 

The formation of the exosphere separates these two eras and 
the regimes of early mass-gain from the nebula and the eon long 
period of particle loss processes. An exosphere exists if there 
is a layer in an atmosphere or generally in a gas sphere where 
the mean free path of the particles is sufficiently large to allow 
them to escape from a planet. Protoplanets fill their Roche-lobe 
because that is where they are in contact with the nebula. For 
a small planet/star mass-ratio the size of the Roche lobe can be 
well approximated by the Hill-radius, Rum = '^(A^pi/[3Mstar])'''', 
where d is the orbital distance. Equating this to the mean free 
path for hydrogen and assuming P = nkT we obtain an exo- 
sphere formation pressure f'exo-form at the Hill-sphere boundary 
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where cth is the hydrogen collision cross section. Assuming an 
exponential decay of the nebula pressure, Po with a time-scale 
Tnebj we calculate the time it takes for the lowest pressure around 
the planet to reach Pexo-form and an exosphere to form within the 
Hill-sphere 
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We take cth = ttOq, with the Bohr-radius, ao as a lower limit 
for collisions under typical nebula-conditions. For typical nebula 



pressures and temperatures at the planet forma tion time, taken 
from the predicted CoRoT-planet population dWuchterl et al. I 
|2006), we arrive at fexo-form ~ 30 - Tmb- With a typical global 
Tneb of ~ 10 Myr representing a compromise between vari- 
ous empirical methods (Tneb ~ 1 - 100 Myr) (see Hillenbrand 
et al. 2005 for a recent update focused on dust IR-emission), 
one obtains average exosphere formation-times, fexo-form ~ 50- 
300 Myr Thus exospheres may form at significant ages after 
the planets' host star arrived at the Zero-Age-Main Sequence 
(ZAMS). 



2. Thermal mass loss 

In order to study the evolution of the maximum possible ther- 
mal hydrogen loss rates from close-in Exosolar Giant Planets 
in a mass range of lO^** g (EGP I) to 10^" g (EGP II) at differ- 
ent orbital distances we estimate the thermal mass loss from the 
evolution of the XUV flux which arrives from the host star at the 
planet's orbit. We use in the present study a scaling law derived 
from XUV flux observations of solar-like G stars with different 
age (Dxuv = 6.13r"Vxuv (iRibas et al. i r2005). where t is the 
age of the star in Gyr, while /xuv = 8.5 x 10""^ W m"^ is the 
flux at 1 AU averaged over the planetary sphere. It is scaled to 
the different orbital distances by using an r^^ dependency. The 
maximum possible thermal mass loss rate Fth can be calculated 
by assuming that an atmosphere contains no efficient IR-cooling 
molecules and by using the energy-limited equation which was 
originally derived for Earth-like planets by Watson et al. (1981). 
Recently, Erkaev et al. (2006) modified the energy-limited equa- 
tion for the application to "Hot Jupiters" 
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where Rp\ and Mpi are the planetary radius and mass, K is a 
factor which is related to the Roche-lobe effect (lErkaev et al. I 
2006) and is 0.5 for the considered planet, rxuv is the dis- 
tance in the thermosphere where the optical thickness rxuv — > 1 
and the main part of the XUV radiation is absorbed. One should 
note that rxuv is much closer to the planetary radius Rpi than 
the distance used by Watson (1981). By using Eq. (3) with rxuv 
of ~ l3Rpi (Yelle 2004) we can calculate the maximum ex- 
pected thermal mass loss rate for HD209458b (Rp\ = 1.43 ^jup, 
Mpi = 0.69 Mjup) at present time. It is ~ 1.7 x lO" g s"' 
and appears to be in good agreement with Vidal-Madjar et al. 
(2003), Yefle (2004) and Tian et al. (2005). The mass loss rate 
for HD209458b at about 200 Myr after its host star arrived at 
the ZAMS is ~ 6.0 x 10'^ g s"'. We note that Lammer et al. 
(2003) applied Watson's assumption for "Hot Jupiters" and over- 
estimated the energy-limited loss rate by an order of magnitude. 
One can see from Table 1 that thermal evaporation may be an 
efficient loss process for close-in gas giants at d < 0.05 AU 
with Mpi < 5 X 10^'^ g (< 0.25Mjup and exosphere formation 
times < 200 Myr If one proceeds to longer exosphere forma- 
tion times the integrated mass loss is gradually decreasing. At 
larger orbital distances, the thermal mass loss is < 1 %. Thermal 
evaporation from close-in lower mass exoplanets may be strong 
enough to remove their hydrogen envelopes from their cores. 
One should also note that exoplanets with orbits < 0.02 AU will 
experience an even higher thermal mass loss due to the higher 
XUV flux and the mass los s enhancement due to the Roche-lobe 
eff^ect (lErkaev et al. II2006I) . 
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d[AU] 


P[d] 


EGPI : L,h [%] 


EGP II: L,h [%] 


0.02 


1 


~ 87 


~ 8.7 


0.05 


4 


~ 14 


~ 1.4 


0.013 


16 


~ 2 


-0.2 



Table 1. Thermal mass loss in % of the initial planetary mass of 
EGP I and EGP II integrated over the history of the stellar system 
for a representative exosphere formation time fexo-form ~ 200 
Myr. 



3. CME-induced nonthermal mass loss 

Taking into account that tidal-locking of short periodic exo- 
planets may result in weaker planetary magnetic moments, as 
compared to fast rotating Jupiter-class planets at larger orbital 
distances (GrieBmeier et al. 2004), Khodachenko et al. (2006) 
found that the encountering CME plasma may compress the 
magnetosphere and force the magnetospheric standoff distance 
down to heights where ionization and ion pick-up of the plane- 
tary neutral atmosphere by the CMEs plasma flow takes place. 
Assuming for the G-type host star of HD209458b the aver- 
age CME occurrence rate as the one observed on the Sun, 
Khodachenko et al. (2006) found that, depending on magnetic 
protection "Hot Jupiters" at 0.045 AU could have lost over their 
lifetime a mass from 2 % up to more than Mjup. For estimating 
the orbital distance at which Jupiter-class exoplanets maintain 
the main part of their initial mass, we calculate the "maximum" 
possible CME-induced mass loss between orbital distances of 
0.015-0.2 AU. 

In the case of our Sun the maximum expected CME plasma 
density hcme at orbital distances < 0.2 AU is estimated from 
the analysis of CME associated brightness enhancements in the 
white-light coronagraph images. By using an analogy between 
the Sun and solar-type stars the maximum CME density depen- 
dence on the orbital distance d from the star can be assumed as 
a power-law (e.g., Khodachenko 2006, and references therein) 
ncMEid) - no(d/doy^-^, which for «o = 5 x 10^...5 x 10'' cm"-' 
and do = 3/?sun gives a good approximation for the values esti- 
mated from the SOHO/LASCO coronograph images. The aver- 
age mass of CMEs, Mcme, is ~ 10'^ g and the average duration 
of CMEs, TcME, at distances (6...10)/?sun is ~ 8 h. The coUi- 
sion rate between C MEs and "Hot Jupiters" can be estimated by 
iKhodachenko et al. (2006.) 
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where /cme is the CME occurrence rate. Acme and 6p\ are the 
CME and planetary angular sizes, is the latitude distribution 
of CME producing regions. Assuming an average stellar CME 
occurrence rate /cme of ~ 3 CMEs per day (Khodacheno et al. 
2006; and references therein), and taking the size, duration and 
latitude distribution of CMEs close to the solar values (Acme - 
njA to 7r/3, fcME = 8 h, = 7^1^), with the angular size 6p\ 
related to the orbital distance between 0.015-0.2 AU of 0.045 - 
0.0034 rad, one obtains a CME colhsion rate/coi of ~ 0.17-0.3 
hits per day. Using these values an average total CME exposure 
time fcME ~ 1.2 X 10'^ s during a period of 5 Gyr is obtained. 

To calculate atmospheric erosion rates we apply the same 
time-dependent numerical algorithm, which is able to solve the 
system of hydrodynamic equations through the transonic point 
for calculation of the upper atmospheric hydrogen density as 
used by Khodachenko et al. (2006) and Penz et al. (2006). 
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Table 2. Maximum CME-induced H^ pick up ion mass loss rates 
per CME collision for "Hot Jupiters" with different substellar 
planetary magnetopause obstacles as a function of orbital dis- 
tance. 



The obtained density profiles are in agreement with the den- 
sity profiles calculated by Yelle (2004) and Tian et al. (2005). 
For studying the expected maximum possible atmospheric ero- 
sion due to CMEs we assume magnetospheres which build ob- 
stacles against the CME plasma flow at substellar distances of 
1.3, 1.5 and 2.0 R-p\. For the CME-induced H^ pick-up loss 
rates caused by the CME plasma flow we take into account the 
ionization processes produced by CMEs during the interaction 
with the neutral atmosphere above the assumed planetary ob- 
stacles. The H^ pick up loss calculations are performed with 
a numerical test particle model which was successfully used 
for the simulation of ion pick-up loss rates from various plane- 
tary atmospheres (L ichtenegger et al. 1995; Erkaev et al. 20051 
IKhodachenko et al. II2006I) . The calculation of the particle fluxes 
is performed by dividing the space above and around the plane- 
tary obstacle onto a number of volume elements AVi. Production 
rates of planetary H^ ions are obtained by calculation of the 
CME plasma flow absorption along streamlines due to charge 
exchange, photo-ionization, and electron impact ionization with 
the extended upper atmosphere. The CME plasma flux Ocme in 
a volume element AVi at a given position with respect to the 
planetary center is determined by 



<l>CME(ri) = ^tL'^r,)e 



(5) 



where the integration is performed from the upstream CME 
plasma flow to the corresponding point s\ at position r\ on the 



streamline. is the unperturbed CME plasma flux, hh is 

the neutral hydrogen density as a function of planetary distance, 
cr^ is the energy dependent cross section of hydrogen ioniza- 
tion processes, and ds is the line element along the streamline. 
The total planetary ion production rate from atmospheric hydro- 
gen atoms in each volume element can be written as the sum 
P\\+ ~ Pn* Pn* /'h+' ^i*-^ Ph+ '■^^ electron impact ionization, 
p'^^ the charge exchange rate, and p^^ the photo-ionization rate. 
The total H^ ion mass loss rate finally becomes 



FcME = ^ /?H+^^i'"H+ 



(6) 



Table 2 shows the maximum expected CME-induced H^ pick 
up mass loss rates due to CME collisions from weakly magne- 
tized "Hot Jupiters" with substellar planetary obstacles at 1.3 
/?pi, 1.5 /?pi, and 2.0 R^x at various orbital distances. Fig. 1 
shows the residual mass after 5 Gyr of the maximum CME 
plasma exposure of weakly magnetized "Hot Jupiters" with 
planetary obstacles at 1.3 /?pi, 1.5 R^x, and 2.0 Rp\ as a function 
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Fig. 1. Residual mass after 5 Gyr exposure of dense CMEs in 
units of Mjup of weakly magnetized "Hot Jupiters" as a func- 
tion of orbital distance compared with planetary obstacles at 
1.3 (solid line), 1.5 (dashed line), and 2.0 (dotted 
line). The horizontal dashed-dotted lines indicate the masses of 
Saturn, Uranus, Neptune, the horizontal dotted lines corresponds 
to the mass of the Earth (lower line) and to the "Super-Earth" 
(Mpi a; lOMEaith) mass domain (upper line), and the black cir- 
cles show the so far observed lower mass close-in exoplanets. 



of orbital distance compared to 5 known low-mass exoplanets 
(Gliese 876d, 55 Cnc e, HD69830b, HD160691d, HD69830c) 
(http://exoplanet.eu/catalog.php). The results obtained from the 
numerical ion pick up test particle model simulations indicate 
that weakly magnetized gas giants at orbital distances < 0.05 
AU or periods < 4 days may lose a mass equivalent to that of 
Jupiter during a 5 Gyr time period. By comparing the results of 
the integrated CME-induced mass loss with the thermal evapo- 
ration loss shown in Table 1, one can see that for magnetically 
unprotected planets nonthermal loss processes are much more 
efficient than thermal evaporation. 

Baraffe et al. (2004) compared the ratio of the mass loss 
timescale ty^ to the thermal timescale fth, which is characterized 
by the Kelvin-Helmotz timescale, found that, when ^/fth be- 
comes < 1 the evolution of the "Hot Jupiter" changes drastically 
which results in rapid expansion of the planetary radius and en- 
hances the mass loss. Depending on the strength of the planetary 
magnetic field, XUV flux and exosphere formation time, close- 
in gas giants at orbital distances < 0.05 AU may experience this 
violent mass loss effect. 

One can also see from Fig. 1 that it is very unlikely that 
the three low mass exoplanets HD69830b, HD160691d, and 
HD69830c are remaining cores of eroded gas giants. They could 
have lost their hydrogen-envelopes but their initial mass was not 
much larger than that of Neptune. The two other low mass exo- 
planets Gliese 876d and 55 Cnc e are located at orbital distances 
where they could have been strongly afi'ected by mass loss. 
Therefore, both exoplanets can be remaining cores of eroded 
weakly magnetized gas giants. 



Gliese 876d and 55 Cnc e, while the three other known 
lower mass close-in exoplanets HD69830b, HD160691d and 
HD69830c belonged most likely since their origin to the 
Neptune-mass domain. The results of our study indicate that 
only a combination of the size detection with CoRoT and 
ground-based follow-up mass determinations together with the- 
oretical mass loss studies over the exoplanets' history can 
bring reliable information on the statistics of remaining cores, 
shrinked gas giants and "unaffected" lower mass exoplanets, like 
Super Earths or Ocean planets. Furthermore, thermal evapora- 
tion from "Hot Jupiters" at orbits < 0.02 AU, together with 
CME-triggere d nonther mal mass loss processes and tidal inter- 
actions ( Patzold and Rauer 2002) can be a reason that so far no 
"Hot Jupiter" is observed at these close distances to its host star. 
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4. Conclusions 



Our study indicates that mass loss of weakly magnetized short 
periodic "Hot Jupiters" can produce the observed masses of 



